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Abstract: Recent ATLAS and CMS measurements show a slight excess in the WW 
cross section measurement. While still consistent with the Standard Model within 1— 
2-0", the excess could be also a first hint of physics beyond the Standard Model. We 
argue that this effect could be attributed to the production of scalar top quarks within 
supersymmetric models. The stops of ~ 200 GeV has the right cross section and 
under some assumptions can significantly contribute to the final state of two leptons and 
missing energy. We scan this region of parameter space to find particle masses preferred by 
the WW cross section measurements. Taking one sample benchmark point we show that it 
can be consistent with low energy observables and Higgs sector measurements and propose 
a method to distinguish supersymmetric signal from the Standard Model contribution. 
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1 Introduction 

The W + W~ diboson production process provides an important test of the electroweak 
(EW) interactions of the Standard Model (SM). Deviations from the SM predictions could 
arise due to new physics contributions, like anomalous triple gauge boson couplings or new 
particles decaying to the same final state as the electroweak gauge bosons. 

The ATLAS and CMS experiments have performed measurements of the WW pair 
production cross section in pp collisions at y/s = 7 TeV and 8 TeV in the fully leptonic 
channel. Using a full dataset at 7 TeV, ATLAS measured the cross section a = 51.9 ± 
2.0 (stat) db 3.9 (syst) ± 2.0 (lumi) pb [1], while quoting the SM prediction at next-to- 
leading (NLO) order of a = 44.7 ± 2.0 pb at = 7 TeV [2]. CMS measurements gave 
a = 52.4±2.0 (stat) ±4.5 (syst) ±1.2 (lumi) pb [3], compared to the SM expectation of a = 
47.0 ±2.0 pb [4]. 1 At yfs = 8 TeV, only CMS has published the results using an integrated 
luminosity of 3.54 fb -1 . It reported a = 69.9 ± 2.8 (stat) ± 5.6 (syst) ± 3.1 (lumi) pb [5] 
compared to the electroweak theory prediction of a = 57.3^'g pb [4]. 

While the above results are far from being conclusive, there is a clear tendency at both 
experiments and center-of-mass energies for a slightly higher measured rate than the SM 
predictions. Interestingly, other EW measurements tend to be in a far better agreement 
with the SM than the WW cross section measurement. This provokes us to speculate that 
the origin of the discrepancy could be attributed to physics beyond the Standard Model 
(BSM). Based on lepton kinematic distributions, ATLAS [1] imposes stringent limits on 
the anomalous WW Z and WW^f couplings. This leaves us with an exciting possibility of 
new particles being produced that contribute to the same final state — two leptons and 
missing transverse energy — as WW pairs. 

Production of supersymmetric (SUSY) particles could significantly affect measurement 
of WW cross section in the fully leptonic final state. It was suggested in ref. [6] that in 
scenarios with charginos as the next-to-lightest supersymmetric particle one could expect 

1 CMS and ATLAS use different methods to calculate the SM cross section, hence slightly different result. 
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an excess in the WW cross section measurement, while avoiding constraints from searches 
in other channels. However, the size of enhancement is limited by the LEP limits [7] on 
the chargino mass. Nevertheless, the chargino contribution can be significant and would 
allow to decrease the tension between the prediction and measurement, provided charginos 
are light and close to the existing bounds, m~± ~ 0(100 GeV). 

The other example of supersymmetric process that could contribute to the WW 
cross section measurement is pair production of top squarks, as we argue in this paper. 
Light stops, motivated by naturalness argument [8-11], are extensively searched for at the 
LHC [12-15]. Cross section is not a limiting factor here — for ~ 200 GeV it easily 
exceeds 10 fb. On the other hand, since stops decay hadronically one has to suppress the 
number of jets in the final state, in order to contribute to the leptonic final state without 
jets. This can be achieved by placing a chargino with a mass only slightly lower than the 
stop mass. The 6-jets produced in the two-body stop decay, t± —> xfb, would be then 
too soft to be reconstructed. The chargino would further decay with on- or off-shell W, 
contributing to the dilepton final state, 

h ^xfb^xl b^xUvb. (1.1) 

The other possibility could be provided by three- or four-body stop decays where kinematics 
also limits pr of 6-jets, however keeping in mind limits from the LHC searches [11, 16, 17]. 
The stop production with a subsequent two-body decay is on the other hand constrained 
by a dedicated ATLAS study [13]. However, because of the applied rriT2 cut, sensitivity 
of this search does not significantly affect a part of parameter space where W becomes 
off-shell. Therefore, in section 3 we fit the signal of the stop pair production, followed by 
the decay chain eq. (1.1), in order to find the minimal supersymmetric standard model 
(MSSM) parameters compatible with the WW cross section measurement. 

The paper is organised as follows. In the next section we briefly discuss the WW 
cross section measurements, the relevant top squark search and simulation procedure. In 
section 3 we perform a scan of the stop-neutralino masses to find a region consistent with 
the WW excess and discuss a method to distinguish SUSY signal from SM processes. 
Finally, we conclude in section 4. 

2 WW and stop searches 

Both ATLAS and CMS have published WW pair production searches. ATLAS measured 
the WW production cross section in pp collisions at yfs = 7 TeV [1] , while CMS published 
results for = 7 TeV [3] and 8 TeV [5] using £ int = 4.92 fb" 1 and 3.54 fb" 1 , respectively. 
As discussed in Introduction, in both cases there was an excess in the observed number 
of events compared to the SM prediction. The experiments were looking at the leptonic 
channel, where the final state consists of two oppositely charged leptons (the same or 
opposite flavour) and missing transverse energy, £ + £~ + E™ ms . In the following we briefly 
recapitulate the ATLAS and CMS searches. 

The main SM backgrounds for pp — > W + W~ —> H, + l~uv process originate from top 
quark production, Drell-Yan processes and other diboson pairs. In order to suppress top 
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quark contribution a jet veto is applied. An event is rejected if there is at least one jet 
with pt > 25(30) GeV in ATLAS (CMS) search. Drell-Yan production is suppressed using 
a cut on the invariant lepton mass, m«, and a projected (relative) EJpf^ defined as 

pmiss _ / ^ iss x sin Afaj if Afa < vr/2 

^T,rel " j E ^ s if A ^ > ^ » l^J 

where A(j>tj is a difference in the azimuthal angle between p™ lss and the nearest lepton 
(jet). After the cuts one obtains relatively clean sample of WW events, with purity of 
~ 70%. The remaining background contribution is estimated using data-driven methods. 2 

Finally, we discuss the search for light stops performed by ATLAS [13], which covers 
a mass region relevant for our study. It targets the same final state as WW analyses, two 
leptons with missing transverse momentum, but using a different set of cuts. Crucially, 
the signal regions in this study require mx2 > 90 GeV. The mx2 variable [18, 19] has a 
sharp kinematic edge at the W boson mass for ti production. For the supersymmetric t\t* 
production the kinematics could significantly differ from that of the top pair production, 
because of an additional contribution to missing energy due to the lightest supersymmetric 
particles (LSP). Therefore, stop production would populate a region of high mx2> where the 
SM backgrounds are suppressed. The situation changes for nearly-off-shell and off-shell W 
in eq. (1.1). In this case, the mT2 cut will also result in suppression of the supersymmetric 
signal and lost of sensitivity. Since ATLAS presented search results for a similar scenario 
with mi x — = 10 GeV we can easily apply those exclusion bounds in our study. 

In order to find a range of stop parameters consistent with experimental searches we 
simulate events using Herwig++ 2.5.2 [20, 21] and process them using fast detector sim- 
ulation Delphes 2.0.3 [22]. We implement selection procedures and cuts for the relevant 
ATLAS and CMS searches discussed above. Furthermore, we validate the implementation 
by comparing efficiencies as reported by the collaborations and we find a good agreement. 
Nevertheless, whenever possible we use the event rates of WW and other SM processes 
given in the ATLAS and CMS notes. The stop signal is scaled to the NLO rate using 
Prospino 2. 1 [23]. With this setup, we perform a scan described in the next section. 



3 Stop contribution 

3.1 Fitting a simplified model 

Given that the stop pair production events followed by the decay chain eq. (1.1) contribute 
to the signal regions of the WW measurements, the following questions should be addressed: 

• Which mass region can fit each experimental result well? 

• Are those mass regions consistent with each other? 

• Are those mass regions consistent with direct stop searches? 

2 At this point the Higgs boson contribution, h — > WW, is not taken into account. 
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• How one can distinguish the stop contribution from genuine WW events? 

Postponing the last question to the next subsection, we address the first three in this 
subsection based on the simplified model approach. 

Our simplified model considers exactly the same process as given by eq. (1.1). As 
discussed in Introduction, the mass difference between the stop and chargino has to be 
small, otherwise the 6-quark from the stop decay would be reconstructed as a high-py 
jet and the event would be rejected by jet veto. We therefore fix the chargino mass by 

— ?n~± = 10 GeV. With this assumption, the model is defined by two parameters: 
and tu^o. As mentioned in the previous section, ATLAS has recently presented the light 
stop search results using exactly the same simplified model. Therefore, one can simply 
apply their exclusion limit to our simplified model parameter space. 

To find out which mass region fits the experimental results, we estimate the x 2 variable 
for each measurement as a function of the stop and neutralino masses: 



Xi(ra t V m x°) = — " ~2— ~ ( 3 - 1 



^I-^sm-^susyKv^x?) 12 



where i specifies the measurement (i = ATLAS 7 [1], CMS7 [3], CMS8 [5]), N® s is the 

(i) (i) 

number of observed events in the signal region, iVo™ and Ag^gy are the predicted contri- 
butions from the Standard Model and SUSY, respectively, and U{ is the quadrature sum 

(i) 

of the systematic and statistical errors. The includes not only the WW contribution 
but also the other SM contributions such as ti and h — > WW* processes. All the factors, 
except for the Ag^ SY , are provided in refs. [1, 3, 5]. 

We estimate ^VguSY( m ti» m x°) m ^ e f°ll° wm g procedure. We generate a grid in the 
(m^m^o) plane with a 10 GeV x 10 GeV step size. In each grid point, 10 5 events of t\t\ 
followed by the decay eq. (1.1) are generated using Herwig++ 2.5.2 [20, 21]. The events 
are processed by Delphes 2.0.3 [22] in order to take detector effects into account. We 
then apply the cuts used in the WW cross section measurement and estimate the efficiency, 
e^m^m^o). The NLO cross section of the stop pair production, a^m^), is calculated 
using Prospino 2. 1 [23]. Finally, the SUSY contribution to the signal region is obtained 

b y Ar susY( m ti' m xS) = £' mt ' °t( m ti) ' t BR (*i ^Xi)] 2 • ^(m^m^o), where C- mt is an 
integrated luminosity. 

Figures 1 (a)-(c) show the y 2 i n the (m^, m^o) plane for the ATLAS7, CMS7 and 
CMS8 measurements, respectively. The area below a black line is excluded by the ATLAS 
direct stop search [13]. In the white top-left region chargino becomes the LSP. Near the 
boundary of the chargino LSP region, the leptons from the X\ ~^ ^ v x\ decay become too 
soft to be detected, leading to AgJj SY — > 0. Therefore in the vicinity of the boundary the 
X 2 approaches to the SM value. 

As can be seen, the best fit regions of the three measurements form a similar arc-shaped 
area, which is roughly symmetric with respect to the dashed green line. The dashed green 
line shows the kinematical threshold of the xf ~^ W^X? decay. In the region above this 
line, the W becomes off-shell and the lepton from the three-body decay, X\ ~^ ^ v x\i 
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Figure 1. The x 2 j e 1- (3-1), distributions in the (m^m^o) plane for each of the measurements, 
ATLAS7, CMS7 and CMS8. In the (d) panel, a combination of all three is shown. Blue areas 
represent the lowest values of x 2 and the region preferred by the experiments. A green dashed line 
indicates kinematical threshold for \\ ~ * W^X? decay. The shaded region below a black line is 
excluded by the ATLAS direct search [13]. A dashed purple line shows a 68% CL region. 



becomes softer as moving away from the line, which in turn requires a smaller stop mass 
to compensate degradation of the efficiency by an enhancement of the cross section. In 
the region below this line, the W from the two-body decay, xf —> Wxl, becomes more 
energetic as moving away from the threshold. This results in degradation of the efficiency, 
because the lepton and neutrino from the boosted W decay are collimated, leading to a 
smaller projected E2p lss . The neutralinos do not contribute much to the E™ 1SS , because 
in the near-threshold region they tend to be back-to-back in the transverse plane and 
their contributions cancel out. In the opposite limit, m^o <C niw, most of the chargino 
momentum is carried by the W and the neutralino becomes soft. 

The dashed purple curves show the 68% CL regions. The regions are somewhat broad 
for ATLAS7 and CMS7. In fact, the SM prediction agrees with the data within \-a 
accuracy for CMS7, therefore adding the stop contribution does not provide a meaningful 



- 5 - 



improvement. On the other hand, the 1-a region for CMS8 is much more localised around 
mjj < 250 GeV and 80 GeV < m^o < 140 GeV. This is because the discrepancy between 
the data and the SM is at about 2-a level and a large stop contribution is required to account 
for the observed excess. Interestingly, for each measurement a large part of the 1-a region 
is consistent with the ATLAS light stop search constraint. Moreover, the preferred regions 
from the three independent measurements are consistent with each other, although two 
of those provide somewhat broad 1-a regions. This agreement is nontrivial since the cuts 
and the center-of-mass energies are different in these measurements. Figure 1(d) shows the 
combined \ 2 °f an three measurements. As can be seen, a significant part of the preferred 
parameter region is consistent with the ATLAS light stop constraint. 

We would also like to comment on the bottom left corner of the parameter space where 
the models are strongly disfavoured by the data. In this region, the contribution from the 
tit* events is too large. This indicates that an analysis similar to the WW cross section 
measurement can also be applied to the light stop search. In fact, the disfavoured region 
spreads to the yet unconstrained area. A dedicated analysis along these lines would be 
able to extend the stop exclusion limits. 

If we are indeed observing the stop contribution, the stop events can fit not only the 
number of observed events after the cuts but also any observed distribution. Therefore 
we compare the data and our light stop model to the distributions provided in ref. [5]. 
Figure 2 shows the distributions of (a) the px of the leading lepton, (b) the px of the 
trailing lepton, (c) the px of the dilepton system and (d) the dilepton invariant mass. We 
choose (to^to^o) = (200,105) GeV as a benchmark point. The histograms show a good 
agreement between the data and the light stop model. The shapes of the SM and the light 
stop contributions are very similar, therefore distinguishing between them would be very 
difficult when using only the provided kinematic distributions. In the next subsection, we 
propose a method to distinguish the stop contribution from the SM. 

The best fit point can be easily realised within the MSSM. With one of the stop states 
heavy, one needs a large splitting between the left and right stops to obtain the Higgs 
boson mass in agreement with experiment [24, 25]. We fix the stop sector by choosing: 
mi R = 195 GeV, = 2000 GeV and At = 2000 GeV. The chargino and neutralino sectors 
are given by: M 1 = 105 GeV, M 2 = 190 GeV, \i = 2500 GeV and tan /3 = 15. Masses of 
other sfermions, Higgs bosons and gluino are fixed by: Msusy = M3 = M40 = 2000 GeV, 
except for the mass of the right bottom squark, mg = 1000 GeV. For such a choice of 
parameters we obtain: = 203.7 GeV, m^o = 104.9 GeV and m^± = 189.5 GeV, in 
the region preferred by the fit. Using FeynHiggs 2.9.4 [26-29] we evaluated the Higgs 
boson mass to be = 125.6 GeV, while the rate in h — > 77 mode turns out to be 
i? 77 = 1.05 • compared to the SM value. Low energy observables have been checked 
withSuperlso 3.3 [30, 31]: BR(B -> A s7 ) = 3.7xl0~ 4 and BR(B S -> /i/x) = 3.45 x 10~ 9 , 
and are consistent with the current experimental values [32, 33]. 

3.2 Stop's smoking gun 

If the excess in the WW cross section measurement is confirmed with a higher significance, 
it will be crucial to confirm that it indeed originates from beyond SM physics. Therefore, 
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Figure 2. Distributions of (a) the leading lepton transverse momentum p™ ax , (b) the trailing 
lepton transverse momentum p™ ln , (c) the dilepton system transverse momentum p^, and (d) 
the dilepton invariant mass m«. The SM, Higgs and stop contributions are shown separately. The 
genuine stop contribution is also depicted for comparison and multiplied by factor 5 for convenience. 
The SM event numbers, data points and errors are taken from ref. [5]. We follow a convention 
proposed in ref. [6] in presenting this plot. 



we discuss here an angular distribution that could help to discriminate between the SM 
contribution and supersymmetric origin. As a working point we choose the benchmark 
scenario discussed in the previous subsection: = 200 GeV, m~± = 190 GeV and 

m^o = 105 GeV. 

Due to different spins and production mechanisms of W bosons and top squarks one 
can expect differences in the polar angle distribution of initially produced particles in the 
hard process center-of-mass frame, as discussed in refs. [34-37]. This indeed is the case 
as can be seen in figure 3(a), where WW production exhibits a strong enhancement in 
the forward direction. In case of stops, the effect is much less pronounced even though 
the forward direction is also preferred. As discussed in ref. [37], such a difference could 
affect angular distributions of the final state particles and provide a strong discrimination 
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Figure 3. (a) The polar angle of the initially produced WW and stop pairs in the center-of-mass 
of the hard process frame, (b) The pseudorapidity difference of the lepton pair, eq. (3.2), for the 
WW and stop signals. 



between different models. 

In order to probe the production distribution more directly, we use the following ob- 
servable [34]: 

cos9* ee = tanh (^^j , ^Vee = Vh ~ % > ( 3 - 2 ) 

where Arju is the difference of the pseudorapidities between the leading and the trailing 
lepton. This variable is the cosine of the polar angle of the leptons with respect to the 
beam axis in the frame where the pseudorapidities of the leptons are equal and opposite. 
Being a function of the difference of pseudorapidities, it is longitudinally boost-invariant. 
Figure 3(b) shows cos 9* u distribution for WW and t\t\ pairs. Much of the difference seen in 
figure 3(a) is now absent, which makes distinction between the two processes significantly 
more difficult. As pointed out in ref. [34], the cos9* u observable requires high boosts of 
initially produced particles. However, both WW and stops have a significant fraction of 
events produced close to the threshold, what partially dilutes the expected difference in 
the final state distribution. 

To improve discriminating power of cos 9$* one should take events with higher center- 
of-mass energy of the hard process. This can be achieved using a variable defined as [38, 39] 

V^min = VE 2 ~ PI + E™ SS , (3.3) 

with E, P z being the total energy and longitudinal momentum of the reconstructed leptons. 
We find that the cut V§ m i n > 150 GeV leads to the highest significance for discriminating 
the WW and stop signals. Figure 4(a) shows cos 9^ distributions after this cut. 

Finally, we discuss the significance of pinning down the alleged stop signal. We follow 
here the approach proposed in ref. [37] and define the following asymmetry: 

A = iV(|cos^| >0.5)-jV(|cos^| <0.5) 
Ntot 
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Figure 4. (a) The pseudorapidity difference of the lepton pair, eq. (3.2), for the WW and stop 
signals after the cut V^min > 150 GeV. (b) The V§ m i n distribution for the WW, tit* and tt events. 
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Table 1. The asymmetry, eq. (3.4), for the WW and tit* pair production without vSmin cu t and 
after applying the V§min > 150 GeV cut. 

where N(. . .) is the number of events fulfilling the respective condition. After applying the 
CMS cuts [5] we obtain the values for A listed in table 1. We compare asymmetry for WW 
and stop production at different center-of-mass energies. Clearly, after application of the 
V§mm > 150 GeV cut we get a better separation of the WW and SUSY contributions. An 
additional cut will decrease the number of events, however as can be seen in figure 4(b) 
more so for gauge bosons than for stops. Therefore, we obtain a cleaner sample with a 
preferable kinematics, so one could expect a better sensitivity. 

Figure 5 shows the expected significance of measuring a difference in the asymmetry 
between SM-only (ie. WW and SM backgrounds) case and SM+t\t*. A clear advantage of 
using VSmm cut is visible. At \/s = 8 TeV, with an integrated luminosity £i nt = 20 fb^ 1 , a 
3-sigma evidence is possible. This can be strengthened when combining data collected by 
both ATLAS and CMS. At yfs = 14 TeV the significance builds up much quicker, providing 
5-sigma discrimination with a few fb" 1 of data. 

4 Conclusions 

In this paper, we have discussed explanation of the excess in the WW cross section mea- 
surement by the production of supersymmetric partners of top quark. The stop production 
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Figure 5. The significance of distinguishing the SM-only and SM-KiiJ case as a function of an 
integrated luminosity using the asymmetry eq. (3.4). The red curve shows the significance with the 
cut Vs m in > 150 GcV, while the black curve without the -\/s m in cut. Different pp center-of-mass 
energies are shown for comparison. 



could provide the right amount of additional signal in the dilepton plus missing transverse 
energy final state, while not producing an excess in the stop searches. The large QCD- 
driven stop cross sections makes it favourable to other possible explanations, like gaugino 
production within the MSSM. The only requirement in the case of stops is suppression of 
jet activity, which can be achieved if the mass difference between stop and chargino is small 
or in the case of three- and four-body stop decays. 

We scan the parameter space of the light stop and the lightest neutralino masses to 
find a region favoured by the present ATLAS and CMS data. The preferred region is 
localised below mj ~ 250 GeV with two branches going down the stop masses. While one 
of them is excluded by the direct stop searches, the other one remains consistent with the 
existing limits. It roughly follows a region where mr = rra^o + mw. In this region, the 
kinematic distributions of the stop signal are very similar to the SM WW pair production 
distributions. We compare the distributions for a chosen benchmark point, obtaining a 
good agreement with the results reported by the collaborations. Finally, it can be easily 
fitted to the Higgs results and the low energy observables. 

If the excess is confirmed with a higher significance in a full 8 TeV data set, it will 
be crucial to establish its true nature. Therefore, we have proposed an observable, based 
on the angular distributions, that could help to distinguish between the SM contribution 
and the genuine stop signal. If the stops are the source of the excess, the full 8 TeV data 
set should provide an evidence for its BSM origin. On the other hand, if the additional 
data do not confirm the excess our results can be translated to the exclusion limits in yet 
unconstrained region of the stop parameter space. 

A final confirmation of the nature of the excess will require more detailed studies. In 
particular, one has to show that the new particles decay to the third generation quarks. 
This task may turn out to be very difficult at the LHC if a mass difference between stop and 
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chargino is very small. In such a case, the final confirmation would require a linear collider 
with a much higher sensitivity to soft objects. It would be a very interesting scenario for 
such a machine, with a few new particles in the kinematical reach, allowing for a high 
precision study of their properties. 
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